Are better conducting molecules more rigid? 



Young-Ho Eom and Hawoong Jeong 
Department of Physics, Korea Advanced Institute of Science and Technology, Daejeon, 305-701, Korea 

Juyeon 

Department of Physics, Pusan National University, Busan 609-735, Korea 

Henri Orland 

Service de Physique Theorique, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France 

(Dated: March. 22, 2006) 

We investigate the electronic origin of the bending stiffness of conducting molecules. It is found 
that the bending stiffness associated with electronic motion, which we refer to as electro-stiffness, 
K e , is governed by the molecular orbital overlap t and the gap width u between HOMO and LUMO 
levels , and behaves as K e ~ t 2 /\/u 2 + t 2 . To study the effect of doping, we analyze the electron 
filling-fraction dependence on K e and show that doped molecules are more flexible. In addition, to 
estimate the contribution of K e to the total stiffness, we consider molecules under a voltage bias, 
and study the length contraction ratio as a function of the voltage. The molecules are shown to be 
contracted or dilated, with K e increasing nonlinearly with the applied bias. 

PACS numbers: 



Conducting polymers have been extensively studied 
over the last decades because of their vast electronic de- 
vice applications. They have 7r-orbital overlapping along 
a conjugated backbone and a gap between the highest 
filled and the lowest unfilled bands, forming a band struc- 
ture similar to inorganic semiconductors. When an elec- 
tron is removed or added, the conductivity is greatly en- 
hanced, allowing these semiconductors to be utilized as 
organic electronic device. On the other hand, conduct- 
ing polymers are much more flexible than semiconduct- 
ing solids, while their electronic functions are very simi- 
lar. This property has led to studies for exploring their 
mechanical deformations such as bending and expan- 
sion/contraction driven by electric triggering. It is known 
that the oxidation causes the expansion or the contrac- 
tion of conducting polymers .depending on whether ions 
are inserted or eliminated Also, chemical doping 

drives electromechanical deformation by delocalizing the 
7r-electrons 0,0]. 

The structural stiffness of a material is determined by 
many factors, such as atomic binding energies, the na- 
ture of molecular bonding, inter-polymer adhesion, to 
mention a few. It is worthwhile to note that the orbital 
stack would be an additional source for mechanical rigid- 
ity, because in order to achieve better electronic conduc- 
tion, structural deformations leading to a loss in orbital 
overlap are less favored. Despite this great diversity of 
causes for the mechanical stiffness of materials, in this 
work, we present the first step to relate electronic prop- 
erties of these materials with their mechanical stiffness. 
This is essential not only to relate the electronic origin of 
stiffness with chemistry-based factors but also to grasp 
the mechanical response to electric signal which differs 
from polymer to polymer. For instance, the extension 



ratio of Polyalkylthiophene (PAT) is larger than that of 
Polypyrrole (PPy) 0. It is interesting to relate it to the 
difference in their electronic band structures. Although 
they have comparable orbital overlap, the gap between 
HOMO and LUMO levels in PAT is about half of that in 
PPy Q. This indicates that not only orbital overlap but 
also band gap comes into play for mechanical stiffness. 

In this work, we investigate the effect of the elec- 
tronic properties of conducting polymers on their con- 
formational stiffness. This enables us to understand 
mechanical responses as well as to estimate the elec- 
tronic contribution to the total structural rigidity. For 
this purpose, we model conducting polymers as a one- 
dimensional chain composed of A-monomer where inter- 
monomer hopping of electrons is allowed via 7r-orbital 
stacking. More specifically, to investigate the bending 
rigidity, we assume the hopping strength to depend on 
the angular configuration of monomers [6]. Furthermore, 
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FIG. 1: Schematic picture of a conducting polymer be- 
tween two electrodes. The one with the anchored polymer 
is grounded and the other is fixed at voltage V. The polymer 
configurations are discretized in a way such that the bond 
angle and the mean position of two adjacent beads define 
the tangential angle <J> n and position of nth monomer, respec- 
tively. 
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to simulate the HOMO-LUMO gap, an alternating on- 
site potential is included in the Hamiltonian. Following 
a semi-classical approach, we trace over the electronic de- 
gree of freedom and obtain the effective potential for the 
angle deformations. We find that the bending stiffness 
associated with electronic properties, which we refer to 
as electro-stiffness, governed by the molecular or- 

bital, t and gap width bet ween HO MO and LUMO level, 
u, and scales as n e ~ t 2 j\Ju 2 + t 2 . Further by analyzing 
the electron-filling fraction dependence on K e , we show 
that doping would make the molecules more flexible. To 
evaluate the specific contribution of n e to the total rigid- 
ity, we consider molecules under a voltage bias and ex- 
amine the length contraction as a function of the bias. 
It turns out that the applied voltage alters the electro- 
stiffness, resulting in the length contraction (dilation) of 
the molecules. 

We consider a conducting molecule in the presence of 
an external electric field. The Hamiltonian for the elec- 
trons responsible for the conducting behavior is taken as 

H e = - y^^„ : » + i(c] t c n+ i+h.c)+y^[(-l)"u-fa: n ]4c K . 

n n 

(i) 

Here the inter-site hopping integral i„, n +i is determined 
by the degree of 7r-orbital overlap, and is maximum when 
the molecule is in a straight form. When the molecule 
is not in the straight conformation, the overlap is de- 
creased, which tends to suppress electron hopping. We 
incorporate this fact by introducing an angle-dependence 
in the hopping integral as t n , n +i = tcos(0„+i — <p n ) = 
tcos9 n - For n — 0, the hopping parameter is spa- 
tially uniform and maximized by parallel orbital arrange- 
ment. In the second term we introduce an alternating 
on-site potential, yielding a gap whose width is deter- 
mined by u. This enables us to investigate semiconduct- 
ing molecule having a gap between HOMO and LUMO 
level. Also, the on-site potential is position dependent 
due to the applied voltage V and v = \e\Va/L is the 
voltage drop across one monomer with electrode spacing 
L(see Fig. 1). Taking the monomer spacing as a, we write 
x n = E"=o cos ^ ~ (1/2) cos 0„. 

The Hamiltonian contains several electronic factors 
that contribute to the structural rigidity: (i)the electron 
hopping favored by a straight configuration; (ii)the band 
gap making the molecules less conducting and tending to 
diminish the effect of the factor (i) ; (iii) the applied volt- 
age bias inducing length contraction or extension. The 
rigidity from electronic origin, K e<n , can be obtained by 

kgT d 2 

« e ,n = y~ Q02 lnTr {c,ct} exp(-/3H e ), (2) 

which can be site-dependent when boundary effects are 
considered. In fact, the rigidity also has contributions 
from molecular bonding and atomic binding potentials, 
which we denote as n m , and the potential governing the 
structural rigidity is written as 



where K n = K e<n + n m . Before proceeding to the numer- 
ical evaluation of Eq. @ , we study the small bending of 
the molecules in the absence of a voltage drop. Expand- 
ing the angle to quadratic order, we define 

^ 0) = -tE^+i+M+^Ef- 1 )"^ ( 3 ) 

n n 

= ^^(ct Cn+1 + h. c ), 

n 

where the position of the nth monomer for small cj>'s 
is denoted by x n . Assuming He to be a small per- 
turbation, we can write —kBTla(e~^ n ^)o ~ (He^)o, 
where (X) = Z^Trfe"^' ' with Z = Tre"^ '. 
It is convenient to work in the Fourier space: c„ = 
iV~ 1/2 £ n e ife "c fc . The Hamiltonian H e 0) can then be 
straightforwardly diagonalized by the canonical trans- 
formation c fc = cosxfca/c,+ + sinxfcafc,- and c k+1T = 
-sinxkak,+ + cos Xkdk- as 

W e °) = J2 

k,a— ± 



where X(k, a) — a^J e 2 (k) + u 2 with e(k) = — 2t cos k and 
tan2xfe = —u/e(k) and a = ±1. Similarly we can write 
Tie 1 ^ in terms of the diagonalizing basis, with off-diagonal 
components. In evaluating (He )o, however, since H^ is 
quadratic in the c's, the only non- vanishing contributions 
can be easily traced and we get 

with the mean number of particles occupying the a = +1 
and —1 bands being determined by 

(A4,a)o = (4 ia a fc ,„> = [e^^-ri + l]" 1 . (5) 

Let us first consider when the system is half-filled so 
that (Afe,+) = 0. From Eq. <@J, it is clear that as 
far as the electronic contribution to the bending stiff- 
ness is concerned, the hopping integral plays a domi- 
nant role: while the denominator in Eq. <@J character- 
izes the band width, the numerator is proportional to 
t 2 . On the other hand, for molecules having comparable 
hopping strengths, those with large band gap would be 
more flexible. The numerical evaluation of Eq. (2J has 
been performed and the resulting stiffness is presented 
in Fig. 2. It is in good agreement with the perturba- 
tive result, Eq. (JJ. While the analytic expression for 
the stiffness can be easily obtained for an infinitely long 
polymer, a finite-sized polymer has boundary effects that 
are characterized by a site-dependent stiffness. As shown 
in Fig. 2, the stiffness is weaker in the polymer bulk than 
in the ends, indicating that when a force is applied, the 
bending of the polymer would be more localized in the 
bulk rather than being uniform all over. 

Even more interesting is the filling factor dependence 
of n e . When hole/particles are introduced in the system, 
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FIG. 2: The electro-stiffness n e vs the hopping integral for 
N = 40 in the absence of a voltage bias. The values of K e 
depend on the monomer position, those near the boundaries, 
e.g. n = 1 being smaller than those at n = N/2. For a given 
n, n e is estimated for u — 20 and u = 60 to show that an 
increase in the gap width suppresses the stiffness. The energy 
scales such as K e ,t and u are measured in units of thermal 
energy at room temperature, fcsT = 0.025 eV, throughout 
this paper. 
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FIG. 3: The stiffness difference A/t e vs the electron filling 
fraction 7 for N = 40 in the absence of a voltage bias: for 
a given n, K e is estimated for t = 60, u = 40 (PPy) and 
t — 48, u = 20 (PAT) to show that the increase of the gap 
width suppresses the stiffness. The inset shows K e for a DNA 
molecule with t = 8, u = 40. Here K e ,t and u are in units of 
thermal energy at room temperature, fcsT = 0.025 eV 



the electro-stiffness becomes weaker than that for a half- 
filled system, as displayed in Fig. 3. We evaluated the 
stiffness difference An e = n e {l) — K e(l/2), where 7 de- 
notes the filling fraction as 7 = N e /N, This is in good 
agreement with the experimental observation that the 
bendingrigidity of Polyurethane film is enhanced by salt 
doping k|. This becomes obvious from Eq. iQJ. When 
the system is half-filled, a = (+1) band would be empty 
so that in Eq. @J, the contribution to the stiffness is 
solely due to the a = (—1) band. On the other hand, 
when holes are doped, the a = (—1) band becomes par- 
tially filled, and the corresponding reduction in (A4,-) 
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FIG. 4: The derivative of the length contraction ratio £ vs the 
voltage bias v for various parameter sets. Here a represents 
the ratio a = l ^ L . The inset shows the length contraction 
ratio I vs v for a — 2. The square symbols and the filled 
squared symbols represent t = 60, u = 40 and t = 48, u = 20, 
respectively. Here K e ,t and u are in units of thermal energy 
at room temperature, fcsT = 0.025 eV 



results in a decrease of K e , The electrostiffness depends 
also on the applied electric field which contributes to the 
energy of the system via a coupling to the charges of the 
electrons: when v -C t and u, we may approximate the 
energy due to the field as Ef — —v J2 n x n(CnC n ) . When a 
positive voltage bias is applied, the molecule increases its 
length, and hence the effective stiffness increases. For a 
negative voltage bias, the molecule tends to contract, re- 
sulting in the reduction of its stiffness. Investigating the 
length deformation as a function of the electric field is 
thus very useful to evaluate the contribution of n e to the 
total rigidity k. As we mentioned earlier, the structural 
rigidity comes not only from those electronic degrees of 
freedom but also from the molecular binding potentials. 
Since for the latter, the deformation is presumably rather 
insensitive to the applied electric potential, the length 
deformation caused by the electric field would directly 
relate the contribution of n e to the total rigidity. To 
see this more clearly, let us define the length contraction 
(extension) ratio £ as 



1 



L 



(6) 



with L c being the contour length of the polymer, and 
(c/> 2 ) being the average angle fluctuation per monomer, 



given by 



= Tr 



{0} f 



-jJH, 



E^f/Tr {0} e- 



-PHr, 



Here, 



disregarding the site-dependence in K ej „ (which as we 
saw is uniform , except for a few boundary sites), and 



writing Ti r , 



E 



j) where k = k t 



and K m is the bending rigidity due to molecular bonding. 
It is clear that (c/> 2 ) oc 1/k and hence, d£/dv, which is 
measurable in experiments, would be simply related to 
the electric-field dependence of K e . Evaluating / as a 
function of the voltage bias for molecules having different 
K m , e.g, a = K m / K e (v — 0, 2, 5), we confirm in Fig. 4 that 
di/dv is indeed not very sensitive to K m . It is also shown 
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that the length contraction ratio increases nonlinearly 
with v and its derivative with respect to v is positive. 
This clearly demonstrates the expected feature that the 
molecules adjust their length to the voltage, allowing for 
their use as electro-mechanical switches. 

In summary, the electronic origin of bending stiff- 
ness was investigated. It was shown that the electro- 
stiffness, K e , is governed by the molecular orbital overlap 
and the gap width between HOMO and LUMO levels: 
molecules with wider band width are more flexible. The 
electro-stiffness can be controlled by molecular doping 
or by applying a voltage bias. Analyzing the electron 
filling- fraction dependence on K e , we showed that doping 
makes molecules more flexible. In addition, we consid- 
ered molecules under a voltage bias to extract the n e 
contribution to the total stiffness. In response to the 
applied voltage, the molecules are contracted or dilated 
with a very nonlinear increase of K e with the applied bias. 

To conclude this study, we mention the value of K e for 
a few molecules. For example, DNA has an extremely 
narrow band width (« 0.01 - 0.04 eV) Q and its elec- 
trostiffness is estimated to be roughly just a few k B T. 
It is well-known that the persistence length of ssDNA 
is about £ p ~ 5nm [§], which is related to the bend- 
ing energy by K — (£ p /a)kBT. Taking the inter-base 



distance a — 3.4, k w 14k bT , showing that the con- 
tribution of K e to the stiffness is significant. For a ds- 
DNA £ p « 50nm and hence n is ten times bigger than 
that for ssDNA Q, while the doubling of n e cannot ac- 
count for the difference. This suggests that among the 
energetic factors which govern dsDNA bending, the elec- 
tronic motions via orbital overlap is not as crucial as the 
electrostatic repulsion between phosphate groups and the 
helical structured base staking @. On the other hand, 
the persistence length of carbon nanotubes (CNT) lies in 
a macroscopic range £ p = 0.1 ~ 1/im ^jj], which shows 
that the bending rigidity of CNT is hundreds times larger 
than that of dsDNA. Noting that the orbital overlap is 
t ~ 2.5 eV, and the band gap is small, u = ~ 0.5eV [IjJ, 
we find K e ~ 10 2 fesT, which shows the important contri- 
bution of electro-stiffness to the total stiffness of CNT. In 
addition, the HOMO and LUMO level of PPy and PAT 
can be simulated by taking u — 1 eV and t — 1.5 eV, and 
u = 0.5 eV and t = 1.2 eV, respectively y|. This leads 
to K e 34fc B T for PPy, and n e w 28k B T for PAT. Al- 
though no direct measurement of the bending rigidity of 
these materials has been made, this goes in the direction 
of showing that the latter is more responsive than the 
former 0|. 
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